A computer controlled pendulum with position readout 
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We have designed, built and operated a physical pendulum which allows one to demonstrate 
experimentally the behaviour of the pendulum under any equation of motion for such a device for 
any initial conditions. All parameters in the equation of motion can be defined by the user. The 
potential of the apparatus reaches from demonstrating simple undamped harmonic oscillations to 
complex chaotic behaviour of the pendulum. The position data of the pendulum as well as derived 
kinematical quantities like velocity and acceleration can be stored for later offline analysis. 



I. INTRODUCTION 

The equation of motion for a simple damped and 
driven pendulum is easy to find, yet its solutions are 
phenomenologically rich. By changing the parameters, 
the behaviour of the pendulum may change from peri- 
odic to chaotic. To demonstrate all possible phenomena 
with a single device it is necessary to apply a well-defined 
and time-dependent external torque at any given angular 
position of the pendulum. 

The idea presented in this paper is to connect the pen- 
dulum both with an angular transducer, which sends an- 
gular position data of the pendulum to a computer, and 
with a torsion spring connected to a stepping motor. The 
angle ■0 between the rotational positions of the axes of 
the stepping motor and the pendulum, respectively, de- 
termines a torque that may depend on the position of the 
pendulum, its angular velocity, and on time. The step- 
ping motor is controlled in a closed control loop formed 
by the angular transducer, a control computer, and the 
motor itself. 

Apart from being used in the described closed control 
loop, the angular transducer also delivers angular posi- 
tion data, with a precision of 8192 steps for one full turn 
of 360 degrees. This data-recording feature allows for 
detailed comparisons with computer simulations^'^'"^'^'— . 

This paper is organised as follows. Section III Al de- 
scribes the mechanical and electronics set up. Software 
concepts are discussed in Sect. Ill Bl In Sections IIIII and 
II VI we collect and discuss several results of measurements 
demonstrating a broad range of applications and the use- 
fulness of the device for teaching the many different as- 
pects of the solutions of the equation of motion for a 
physical pendulum. 



II. THE APPARATUS 



A. Mechanical Set up and Electronics 




The mechanical set up consists of an aluminium base 
structure holding the pendulum, the angular transducer, 
the stepping motor, and the torsion spring. Figure [1] 
shows a photograph of the pendulum; a schematic repre- 



FIG. 1: Photograph of the pendulum. 



sentation is shown in Fig. ^a,). 

The body of the pendulum is made from stainless steel 
and consists of a rod with a length of 0.20 m and a set of 
steel disks to adjust its mass and moment of inertia. By 
changing the number of steel disks the total mass of the 
moving part can be adjusted from 0.110kg to 0.288kg 
and the corresponding moment of inertia varies between 
1.56 X 10"'^ kg m^ and 6.61 x 10~^kg m^. The steel sus- 
pension axis divides the pendulum rod in a length ratio 
of one to three. It is mounted with ball bearings in the 
aluminium base structure. The axis is connected to the 
angular transducer. The transducer works optically with 
a minimum of internal friction. It provides a precision of 
8192 steps for a full turn of 360 degrees corresponding to 
0.77 mrad per step. The resulting internal friction of the 
pendulum is caused mainly by the angular transducer 
and by imperfections of the ball bearings holding the 
steel axis in the base structure. The position data of the 
transducer are sent to a digital I/O-card (DAQCard"^^- 
DIO-24, National Instruments) placed in the PCI port of 
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FIG. 2: Panel (a) shows the mechanical setup of the pendu- 
lum, (b) the definition of the absolute position <^ of the pen- 
dulum and of the angle between the axes of the pendulum 
and the stepping motor. 



the control computer. A program converts the grey-code 
of the optical angle encoder into a numerical value for 
the angle ^{t), representing the actual angular position 
of the pendulum at time t, as shown in Fig. [HJb). The 
angular position has a periodicity of 27r. 



The stepping motor allows driving and holding the mo- 
tor axis at a well-defined angular position. The employed 
hybrid motor divides one full rotation into 400 digital 
steps corresponding to 15.71 mrad per step. The signals 
for controlling the motor are generated on a drive board 
connected to the digital I/O-card of the control com- 
puter. The torsion spring is made of a thin steel band and 
was manufactured by a watch-maker. Its spring constant 
has been determined to 0.251 Nm/rad. 



B. Software Concept 

1. The equation of motion 

According to Newton's second law, the sum of all act- 
ing torques is equal to the product of the moment of 
inertia / and the angular acceleration (f of the pendu- 
lum: 

Iip = T* {ip,ip,t) - Tr^niax sin - bip-c sgn((^). (1) 

The sum of all acting torques receives contributions from 
the torque T* generated by the torsion spring via the 
stepping motor, the restoring torque — T^jinax sin of the 
gravitational field, a linear friction term —b (p (caused by 
the air and the control system) and a constant friction 
term — c sgn((^) (caused by the internal friction of the 
ball bearings and the rotation transducer). The value of 
the function sgn((^) is defined as 1 for ip > 0, as —1 for 
ip < 0, and as for 99 = 0. 

The user is able to vary the torque by controlling the 
stepping motor. The computer control of the motor has 
been constructed such that the torque generated by the 
motor has the following functional form 

T* p,t)^ T* sin [ut) + T*^,, sin ^ + 6> + c* sgn{p) . 

(2) 

Here the angular frequency co and the set of four param- 
eters, TQ,r*^3^x,6* and c* are adjustable by the user. 
Inserting the parametrisation of the torque T*, Eq. 
into Eq. JT]) yields the equation of motion 

Ip = -Tr^cB sin ip + Tq* sin(cjt) - bcsp ~ Ccssgn{p>) (3) 

where 



bcS = b - b 
Ccff = c - c 



(4) 



The physical meaning of the four terms on the r.h.s. of 
Eq. ([3]) is as follows. The first term represents the ef- 
fective gravitational acceleration, the second term is the 
harmonic driving term, and the last two terms represent 
the effective friction. The linear or viscous friction term 
(X p is governed by the parameter bcs the (piecewise) 
constant friction term oc sgn(i^) by the parameter Ccff. 
Changing the parameters marked with * the software 
control of the stepping motor allows one to set the ef- 
fective parameters for gravitational acceleration, driving 
and damping to any value required. This key feature of 
the system even allows one to compensate for the friction 
terms completely or to set the effective acceleration of the 
pendulum in the earth's gravitational field to any value 
including zero. In the following we shall, however, always 
suppose that the effective friction parameters satisfy 



bcS > 0, 

Ceft > 0. 



(5) 
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2. The control system 



III. RESULTS FOR A NON-DRIVEN 
PENDULUM 



The software which has been written to achieve this 
flexibihty basicaUy consists of a closed control loop, 
which sets the correct position of the stepping motor at 
any position and velocity of the pendulum at any time. A 
graphical user interface (GUI) allows convenient control 
of the parameters of the equation of motion, Eq. ([3]), for 
setting up the experiment before or even during motion. 

The angular distortion tp of the spring depends linearly 
on the torque T*, 



Do 



(6) 



with the inverse of the spring constant Dq as proportion- 
ality factor. This relation defines the angle ip between 
the axes of the stepping motor and the pendulum which 
is required to obtain the desired torque. 

The actual position of the stepping motor is deter- 
mined by the closed control loop. In each run through the 
loop the actual positions of the pendulum and the step- 
ping motor are compared and, if necessary, the stepping 
motor is moved one step towards the required position. 
The speed of the closed control loop is a linear function of 
the last measured deviation between the actual and the 
required position of the stepping motor as long as the 
maximum stepping rate of the stepping motor is not ex- 
ceeded. The maximum stepping rate is 1300 Hz. As indi- 
cated in Eq. ^ the closed control system also induces lin- 
ear and constant friction terms. Their strengths depend 
on the parameters b* and c* , respectively. Through the 
closed control loop these parameters can be adjusted by 
the user via the GUI. A special method of measurement, 
which is discussed in Sect. IIIIl allows one to find out the 
actual values of these constants b* and c* . The angular 
velocity of the stepping motor is determined by the an- 
gular velocity of the pendulum ip and the corresponding 
angular velocity of the motor axis as required to produce 
the correct angle ip for the desired time-dependent torque 
T*{(f, (fi, t), related to ijj by Eq. ©. The parameters have 
to be chosen such that the maximum stepping rate of the 
stepping motor is not exceeded. A critical parameter for 
the set up of an experiment is the limited maximum hold- 
ing torque of the stepping motor (0.494 Nm) to which 
the mass and the maximum acceleration of the pendu- 
lum have to be adjusted. The graphical user-interface 
has been programmed in C-l — h using the QT- and Qwt- 
libraries for the generation of graphical elements. The 
parameters of the closed control loop, as well as the pa- 
rameters of the equation of motion, can be selected with 
the help of various sliders. Several options for measure- 
ments can be selected by the user. All recorded position 
data are saved and can be plotted in different represen- 
tations. 



Without external driving torque, Tq* = 0, the oscilla- 
tion of the pendulum exhibits specific different charac- 
teristics depending on the selected damping terms. 



A. Small elongation and linear friction 

In this subsection we suppose that the parameters of the 
pendulum are chosen such that 

Tr^c« > 0, 

bcS > 0, 

Ceff = 0. (7) 

We also suppose that the amplitude of the oscillation is 
much smaller than 7r/2. Then the nonlinear equation of 
motion, Eq. ([3]), can be approximated by using simp ^ ip 
which leads to 



Here we define 



ip H — ip + ujQip = 0. 

T 



21 



(8) 



boS ' 



(9) 



The solutions of Eq. ^ are discussed in many 
textbooks^. We have to distinguish three cases. 

(i) Underdamping (oscillatory motion) 
This happens for 

1 



Wo > ^2 ■ 



(10) 



The general solution reads for this case 

(p{t) — (fio exp{—t / t) cos{il!t + S) (11) 

where 



(12) 



The free parameters of the solution, Eq. (fTT|) . are 
ipo and S. The realization of such a pendulum mo- 
tion is shown in Fig. EJJa). 

(ii) Overdamping 
This happens for 



Wo < ^. 



(13) 
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Here the general solution of Eq. ([8]) reads 

fit) = <fi+ exp 
+ ip- exp 




(14) 



with (p^ and (p^ free constants. An experimental 
realization of this case with our apparatus is shown 
in Fig. Hb). 

(iii) Critical damping 
This happens for 



^0 = ^■ 



(15) 



The general solution of Eq. ^ for this case reads 

(16) 

where ipi and ip2 are free constants. 



ip{t) = ( (ySl + (^2^ j exp f 



B. Small elongation and piecewise constant friction 

Here we suppose that the parameters of the pendulum 
have been chosen such that 

Tr^cS > 0, 

bcB = 0, 

ccff > 0. (17) 

For small oscillations we thus obtain, from Eq. ([3]), the 
equation of motion 

pit) + ujI p{t) + ujIxo sgn{pit)) = (18) 

where is as in Eq. ([9]) and 

(19) 



Coff 

X0 = -2- 



For initial conditions which allow for an oscillatory so- 
lution of Eq. (fTS]) this solutions can be constructed in 
the following way. Writing the elongation and its time 
derivative at the initial time i, as 



(pi = p{U) 

a solution of Eq. is given by 



(20) 



= "Xo sgn((^i) + [pi + xo sgn(95i)] cos [uj^it - U)] 



— Pi sin [wo(t - ti)\ . 

Wo 



(21) 
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FIG. 3: Different forms of motion of the non-driven pendulum 
at small elongations. Panels (a) and (b) show results where 
only the linear friction term is present, see Eq. ((Tjl. Panel (a) 
realizes the underdamped, panel (b) the overdamped case, see 
Eqs. and (|14|) . respectively. Panel (c) shows the case of 
small elongations with a piecewise constant friction term, see 
Eq. (|34p . There the amplitude of the oscillation decreases 
linearly with time. 



Inserting the second time derivative of this expression 
into Eq. (fT5|) one finds that 



p{t) + ujI p{t) + Lulxo sgn((/)(i)) 
= WqXo [sgn((^(t)) - sgn((y9,)]. 



(22) 



Hence, Eq. pT|) solves the equation of motion p8)) only 
as long as the angular velocity does not change its sign. 
We illustrate in the following how Eq. (PT|) can be used 
to iteratively construct a damped oscillatory solution of 
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Eq. nil). 

We choose initial conditions where for <i = the pen- 
dulum starts at an elongation 



with infinitesimally negative velocity 



(23) 



(24) 



Here and in the following e > denotes an arbitrarily 
small positive number. For this choice Eq. (PT|) implies 
that 



= Xo + {fo - Xo) cos(woi), 



sgn{(p{t)) = -1 for < t < 



UJo 



(25) 



It is easy to check that Eq. ((25)) gives indeed a solution 
of Eq. dlHl) for 



<t < tt/loo, 



(26) 



that is, until the pendulum reaches its first turning point. 
If (po > 2xo this turning point is located at the negative 
angle 



(1) o 



(27) 



Now suppose that the initial elongation was cpo > 3xo, 



I.e., 



> Xo- 



(28) 



In this case we can continue to construct the solution of 
Eq. (fTSl) . for a further half oscillation period, by setting 
in Eq. ^ 

( (1) 



For times 



TT 27r 
— <t< — 

ujQ Wo 



(29) 



(30) 



the solution is given by 



"^(i) = -Xo - {v'i^ ~ Xo) cos([Joi - 71"), 

^{t) = {'p'a^ - Xo) t^osin(cjoi - tt), 
sgn((^(i)) = +1. 



(31) 



Clearly, we can continue this game, provided that we 
had chosen ip^ > bxo initially. For the elongations at 



the turning-point times, t = mr /loq (n = 0, 1, 2, . . . ), we 
obtain 



^( — ) =(-l)"[^o-2nxo]. 



(32) 



That is, we find a linear decrease of the maximal elonga- 
tions with n. The above way of constructing a solution of 
Eq. ceases to work, though, once the absolute value 
of the elongation at a turning point is smaller than xo- 
We define, therefore, Umax as the positive integer with 
the following property 



-, ^ V'o - Xo ^ 
2X0 



(33) 



The solution of Eq. (jlSp with the initial condition (j23p . 
(I24l) then reads 



ip{t) = (-l)"|xo + [po - (2n + l)xo] cos(wo< - ri7r)| 
(n + 1)tt 



mr 

for — <t < 

UJo 



-, n = 0,l...(n„,ax- 1). (34) 



An experimental example of such a motion of the pendu- 
lum is shown in Fig.ISjc) where the linear decrease of the 
maximal elongations, as given theoretically in Eq. (|32|). 
is clearly visible. 

But what happens with the solution ([M)) for t > 
rimaxTi'/'^o? If <i5('T-max7''/'^o) — 0, the equilibrium has 
been reached precisely at time t — rimaxTr/ajo. It can 
be shown that in all other cases the solution of Eq. (fTSl) 
given in Eq. (|M1) cannot be extended for t > nmaxTr/wo- 

This simply means that for very small amplitudes ip 
and velocities if the approximation of the friction term 
by a function Cefisgn((p) ceases to make sense. In reality, 
the discontinuous function sgn{(p) must then be replaced 
by a continuous function, for instance one behaving as ip 
for 1^ ^ 0. For very small amplitudes and velocities the 
motion of the pendulum will be described by one of the 
cases discussed in Subsection IIII Al That is, the pendu- 
lum will return smoothly to its equilibrium position. 

From data like those shown in Fig. [31 the numerical 
values of the effective friction parameters 6cff and Ccft in 
Eq. ^ can be determined. First, the constant friction 
term is varied until the amplitude approaches an expo- 
nential decrease. This allows to read off from the data 
bcs/I and, since / is known, bcs- Then the linear friction 
term is varied until one reaches the undamped oscillation 
case, that is, absence of the linear friction term. Then 
the constant friction term can be tuned and determined 
from the data. This method of measurement is supported 
by theory only as long as the amplitude of the oscillation 
is small enough. Nevertheless the results are quite satis- 
fying. 



C. Large elongation and no friction 

Now we turn to the case of zero effective friction terms 
which can be easily realized by tuning the torque T*, 
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Eq. ([2]), appropriately. We abandon the approximation 
of small elongations of the pendulum and consider the 
exact equation of motion, Eq. ([3]), with T^^ofr > 0, Tq* = 0, 
^cff — and Ccff — 0: 



where again 



ip{t) +ul sin ip{t) = 



^0 = ^r,cff / i- 



(35) 



(36) 



To solve ([35|) we take as usual^ the route over the energy 
theorem. This gives 



1 



ip^ — cos (f = rj = const. 



where, clearly, we must have 

We shall discuss here only the case 



and set 



ujI<i]<ujo 



rj = -t^o cos(^o, 

< (y90 < TI"- 



(37) 

(38) 
(39) 

(40) 



Denoting again the initial conditions at t — ti as 
<p(^i) = fi, 

^iU) = (41) 

the solution of Eq. (|55|) in determined, via Eq. ([57)l . by 
the integral equation 



dip' 



v^ Jsin^^-sin^^ 



sgn(vjO 2loo {t - U). (42) 



Here 



2 



sin^^ = sin^^ + ^(viO^ 



(43) 



and the initial conditions must be such that Eq. has 
a solution with < ip^ < 'k- The integral in Eq. ([1^ 
is of the elhptic typei. The pendulum oscillates with a 
certain period T between pq and —pa- Contrary to the 
case of small oscillations, the period T depends now on 
the amplitude pq. From Eq. (|^^ we find easily 



9 rVo 



dp' 



'^0 Jo , /sin^ ^ - sin^ f 



(44) 



This can be related to the complete elliptic integral of 
the first kind^ 



1 = — K sm — 

ijJn \ 2 / 



4 
4 



da 



'^o-'o ,/l -sin^^sin'a 



(45) 
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FIG. 4: Dependence of the period of oscillation T, Eq. (|45|l , on 
the maximum elongation po. The data points are compared 
to the calculation which is represented by the solid line. 



Figure [4] shows the calculated dependence of this oscil- 
lation period T on the maximum elongation p^ and the 
verification of this curve by five measurements. 

The period of oscillation also depends on the value 
of the effective gravitational restoring torque Tr^off, see 
Eq. (pS)) . This can be varied by changing the mass of 
the pendulum, which is quite easy to do, or the effective 
acceleration due to gravity, which seems to be impossible 
without leaving the earth or at least its surface. For our 
pendulum we can vary the effective gravitational torque 
7r,eff by varying the parameter T*^^^^^, see Eqs. ©-(H]). 
At each position of the pendulum the torsion spring adds 
a well-defined torque. In this way it is possible to demon- 
strate an oscillation like that which would occur on the 
moon or without any gravitational field. In the latter 
case we have loq = and the pendulum just moves with 
constant angular velocity. This case corresponds to the 
situation of a pendulum in empty space or in a freely 
falling laboratory in a gravitational field. That is, we 
can simulate the conditions prevailing in a space labora- 
tory like the ISS. 



IV. RESULTS FOR A DRIVEN PENDULUM 

A. Resonance phenomena 

Here we consider an external driving of the pendulum 
in a periodic fashion but with small enough amplitude of 
the driving torque such that the pendulum never reaches 
the top position. This leads in general, after the tran- 
sients have disappeared, to an oscillation with frequen- 
cies related to that of the periodic driving torque. In 
this regime we can study resonance phenomena with our 
device. Tuning the constant friction term in Eq. ([3]) to 
zero, Ccff = 0, and supposing Tr^eS > and bcs > our 
pendulum's motion is governed by the equation 



ip H — p> + Wq sm p 

T 



Tr 



'j-sm{ujt). (46) 
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FIG. 5: Amplitude (po and phase shift 5 of a pendulum driven 
by a periodic torque of the frequency The curves are the 
analytical results of the linearised equation of motion, see 
Eq. (|49p . The curves differ from the data points for the the fre- 
quency range of large amplitudes due to the non-linear terms 
in Eq. (ge)). 




tis) 

FIG. 6: Elongation ip{t) of the externally driven pendulum as 
function of the time t. For time < f < 15 s the effect of the 
transients is clearly visible. 



is, the terms cx exp(— t/r) in Eq. (I48p . have decayed the 
maximal elongation of the pendulum, ipo, and the phase 
shift, S, can be measured as functions of the driving fre- 
quency V — Lo I (27r). The measured points shown in Fig.O 
are compared to the results of the calculation for small 
elongations in Eq. P5|) . The data points in Fig. [5] follow 
nicely the theoretical curves except at the peak of the 
amplitude in the resonance region, that is for v 1 Hz. 
There the amplitude 939 gets large and the approximation 
sin Lp ^ becomes bad. 

In Fig. [5] we show results of measurements for ipit) 
for the driven pendulum starting at t = at = with 
some initial velocity. The effect of the transients is clearly 
visible for < t < 15 s. After this time the pendulum 
oscillates with the frequency of the external torque. 



Here Wq and r are as in Eq. 



We shall discuss here 



only the theory for small elongations. Approximating 
sin (/? by iy9 in Eq. (|46p we get 



ip- 



-ip- 



ulip 



-sin(wi). 



(47) 



Again, as in Sect. IIII Ai we have to distinguish the un- 
derdamped, the overdamped, and the critically damped 
cases. As an example we give the well-known solution for 
the underdamped case 



(p{t) — ipo exp(— t/r) coa{ujt + S) 
+ ipo sm{ujt + S) 



(48) 



where uj is given by Eq. (|12p and 



Tr 



tan (5 — 



2lu 



,2\2 



T2" 



^1/2 



(49) 



The constants <pq and 6 are arbitrary and are needed to 
fulfil general initial conditions. After the transients, that 



B. Chaos 

Perhaps the most interesting phenomenon of a driven 
pendulum is chaotic motion^ . One peculiarity of chaotic 
motion is that the oscillation of the pendulum will never 
become periodic, in spite of the fact that it is driven by 
a periodic force. As long as the amplitude of the driv- 
ing force is small enough, the pendulum behaves like a 
harmonic oscillator as discussed in Sect. IIV AI By in- 
creasing the amplitude of the driving force the behaviour 
of the pendulum becomes non-harmonic. In this state 
it is possible to observe both periodic and chaotic be- 
haviour depending on the values of the other parameters 
in the equation of motion. 

Figure [7] shows a nonlinear but periodic oscillation, 
while the motion of the pendulum in Fig. [8] is chaotic. A 
chaotic oscillation can only be reached if the pendulum 
turns over the top position at least once. A critical situ- 
ation is reached when the pendulum just reaches the top 
with zero velocity there. If the pendulum stops infinites- 
imally before the top it will fall back to the side where 
it came from. If the pendulum stops infinitesimally after 
reaching the top it will continue its motion on the other 
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FIG. 7: Elongation ip as function of the time t for a nonlinear 
but periodic oscillation of the pendulum. 



FIG. 9: Phase-space representation of a chaotic motion of the 
pendulum. 
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FIG. 8: Elongation as function of the time t for a chaotic 
motion of the pendulum. 
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FIG. 10: Chaotic motion of the pendulum in phase space 
represented by a Poincare map. 



side. Then, the dependence of the motion of the pendu- 
lum on the initial conditions is so strong, that predicting 
the long term behaviour is practically impossible. This 
is the regime of deterministic chaos^. 

To analyse the nonlinear behaviour of a driven pen- 
dulum it is useful to study its motion in phase space. 
This is here the two-dimensional space of position </? and 
canonically conjugate momentum, in our case equal to 

the angular velocity. Our device allows us to plot the 
measured values of the angular velocity of the pendulum 
vs. its position for each time step. Figure [5] shows the 
phase-space diagram of a chaotic oscillation. The differ- 
ences in the phase-space representations of periodic and 
chaotic oscillations are more obvious if the continuous 
curve is turned into a discrete one. This can be done by 
mapping only the points (<^(i)j V^^S) fo'" times t where the 
phase of the driving torque has some fixed value. This 
is equivalent to a stroboscopic view of the motion of the 
pendulum in phase space at the frequency of the driving 
force. The result is called a "Poincare map." It is ob- 
vious, that the number of points of the Poincare map is 
limited if the motion of the pendulum is periodic. The- 
ory predicts that the Poincare map of a chaotic oscillation 
should show fractal characteristics^. Figure [TUl shows a 



Poincare map with 4000 data points from a chaotic mo- 
tion of the pendulum. Our experimental results are in 
complete agreement with computer simulations made by 
one of us (T. G.) following Baker and Gollub^. 



V. CONCLUSION 

We have constructed and operated a versatile device 
for demonstrating the rich range of motional phenomena 
possible with a damped and driven physical pendulum. 
The degree of damping, restoring, and driving by an ex- 
ternal periodic torque, can all be determined by a com- 
puter controlled torsion spring. A computer readout of 
the position data allows for detailed ofHine study of the 
motion. The device can serve as a teaching tool at schools 
and universities to demonstrate all aspects of oscillatory 
motion ranging from simple harmonic motion to chaotic 
behaviour. 
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